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Some nonlinear aspects of a cluster phenomenon in nuclei are considered using of cubic 
Nonlinear Schrodinger Equation and Korteveg de Vries Equation. We discuss the fol- 
lowing possible nonlinear effects: i) the decribing clusters as solitons; ii) an anomalous 
large angle scattering of a-particles by light and intermediate nuclei; iii) stable vortical 
objects; iv) and dynamical clusterization in the presence of instability 

1. Motivation 

Clusters is a very general phenomenon. They exist not only in nuclei and atoms, 
but also in subnuclear and macro physics (from small drops to atmosphere, stars 
and galaxes) (see matherials and the summary talk of Prof. W. Greiner at CLUS- 
TERS'93 Conference (Santorini,Greece)0 

There exist very different theoretical methods, developed in these fields. How- 
ever there are only few basic physical ideas, and most of methods makes deal with 
nonlinear partial differential equations. 

The aim of the present paper is to show that very different effects of nuclear 
cluster physics could be explained by using of cubic Nonlinear Schrodinger Equation 
and Korteveg de Vries Equation. 

The papers is organized as follows. 

In Sec. 2 we consider general aspects of nonlinear approach for nuclei. An anoma- 
lous large angle scattering of a-particles by light and intermediate nuclei are dis- 
cussed. 

A way to descrive clusters as stable nonlinear localized waves (solitons) is dis- 
cussed in Sec. 3. Section 4 is devoted to stable vortical exitations. Dynamical 
clusterisation and break-up of nuclear system are considered from the point of view 
of nonlinear dynamics in Sec. 5 Last section is a short summary. 
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2. Nonlinear Approach for Nuclei 

The general philosophy giving birth to the NoSE is a rather simple one. Previous 
studiesa into the anomalous large angle scattering (ALAS) of a-particles by light 
and intermediate (a-clustcr) nuclei that the rise in da/dfl (elastic) for backward 
angles - together with the measured resonances, Ex, and decay widths, T^, of 
the dinuclear (a-cluster) system - may be described by (the same parameter of) 
a crude effective surface potential (ESP) arrived at by supplementing the optical 
model interaction, Vqm, by a repulsive (hard or soft) core, Vr'- 

Vesp — VoM + Vr 

VoM facilitates the correct evaluation of da/dVl for forward angles while K gives 
the additional potential contribution allowing for an appropriate description of the 
measured E^, and of the backward rise in da/dfl as observed for 0cm > 90°. 

Similar quasi-molecular features as in these a-scattering experiments have also 
been measured in elastic collisions between heavy ions. Knowing that semi-classical 
and in particular hydrodynamical approaches work nicely when applied to heavy 
ion physics, it appears very suggestive to try to model the elastic scattering of two 
heavy ions in terms of colliding liquid drops with diffuse surfaces. 

Assuming that the "membranes "/surfaces of the two drops remain impenetrable, 
one thus obtain automatically a compression of the densities in the surface regions 
of the two touching drops. The physically motivated boundary condition is then 
that this compession — giving rise to a repulsive spring-force proportional to Vr — 
vanishes for separations 

r > ro = i?i + i?2 + (ai + a2)/2 

between the two ions/drops, i.e. Vr 0. Ri and denote radius and diffuseness 
of the i-th ion; i — 1,2. Thus the dynamical interaction between the two nuclei is 
for r > To characterized by their unperturbed/uncompressed asymptotic densities 
poi{Ri,ai) — corresponding to the measured (charge) densities — and for r < rg 
by the perturbed/compressed densities Pi{r; Ri, ai). 

The central equation emerging from such considerations together with the use 
of the Euler equation and some algebra in liaison with the neglect of higher order 
derivatives is the nonlinear Schroedinger equation (NoSE)0 

_ V2* + V'J'-C(l--^)* = £;*, Po = Poi+Po2 (1) 

2m Po 

where we use the conventional notation with C — K/9 for the nuclear compress- 
ibility, p and Po denote the perturbated and unperturbated densities respectively. 
(Usually K is taken to apply to the bulk of nuclear matter while we consider only 
compressions in the surface regions.) 
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It turned out that Eq. (1) is a narticular solution of a general equation suggested 
by Bialynsky-Birula and MysielskyQ. Using instead of cubic term e log(a 1 1/' | )V' they 
formulated the nonlinear quantum mechanics without any contradiction in classical 
quantum mechanics. By the way, Gaehler, Klein and Zeilingem attempted to find 
the constant e according to the predictions made in the worko. But it was made 
without any success. Our approach supposes the constant C = K/9 or e ( here 
both constants are the same) to be the compressibility constant. It was confirmed 
from the values of the constant K, obtained from different experiments with the 
help of this approach (K « 270MeV) and also by the recent work of J.Braechei^. 
J.Braecher considered the nonlinear term in the NoSE from the point of view of the 
information theory. According to the work of Shannon and Weaver the information 
/ acquired upon measurements of the state | -0 > is proportional to the logarithm 
of the probability I = — log2(a \k \ ) bits, where a =| -00 |~^- 

The certain amount e of energy per bit is expended, transferred, stored, associ- 
ated with the information encoded. That's why it is necessary to add the term el 
to the standard Hamiltonian: 



iJ = i/o + e/ = ^_V2 + C/-elog2(a| VP) (2) 

3 = 1 ™ 

Now it is possible to show that this constant e is very close to our compressibility 
constant. The total energy as the expectation value of the Hamiltonian (2) is: 

E=<i>\H >^< V I iJo I > +e < V' I ^^0 I V >= 

^Eoe<ij\ log2(a | V H | V > (3) 
Defining the entropy for the subsystem: 

5 = <^A|ln(a|0n |0> (4) 

Where K - the Boltzmann constant. Now we obtain the free energy Eq, the tem- 
perature T and constant e : 

Eo = E-TS, T=(^||^^, e = KT\n{2) (5) 

The value e is very close to the constant C. If we suppose the value KT for the 
average kinetic energy of nucleons w SOMeV we obtain the reasonable value of 
compressibility constant. 

Besides NoSE there is another interesting equation, which was applied many 
times for the prediction of new forms of the motion in nucleiEl. This is well-known 
Korteveg de Vrios En uation [KdY) equation. KdV equation was applied many times 
in nuclear physicaj|j-E£lOll3^t3. 

In the works the connection was established between the solution of the KdVE 
and the nuclear potential in Schroedinger equation. At the certain circumstances 
we can see the evolution of the nuclear potential with time. 
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3. Nonlinear nuclear vortex 

In heavy ion coUisions a locahzed stable or long lived excitation might develop. 
The existence of a fairly well-defined nuclear surface makes it possible to discuss 
localized excitations in the nuclear surface region. 

G.N. Fowler, S.Raha and R.M.WeinerM attempted to predict the localized stable 
surface excitations in peripheral heavy ion collisions like in Rossby waves. 

In the workEl A.Sandulesky and W.Greiner obtained a new coexistence model 
consisting of the usual shell model and a cluster model, describing a soliton moving 
on the nuclear surface. They got the excellent description of experimental spectro- 
scopic factors for cluster decays. 

In the framework of semiclassical nonlinear nuclear hydrodynamicsEl we hav 
presented a type of vortical excitations new for nuclear physics. In the paper 
we began to consider a pure vortical motion in nuclear systems excluding the usual 
approximation of smallness of the excitation amplitude and the additional assump- 
tions about the shape of_the nuclear system. We considered the two-dimensional 
analog of a nuclear disk 113, a plane nuclear vortex, a new type of a pure vortical 
state of incompressible inviscid nuclear matter. It is a finite area of constant vortic- 
ity in a plane within a uniformly rotating contour. These states can be considered 
as the generalization of the elliptic Kirchhoff vortex tB. Recently we generalized 
this type of solution to the similar three-dimensional picture and consider a new 
possible class of nonlinear Euler type equations on a nuclear surfacal3. We have 
found a class of uniformly rotating solutions to nonlinear Euler type equations on a 
sphere. The equations of motion describing a localized vortex on a spherical nuclear 
surface - bounded region of constant vorticity, surrounded by irrotational flux - are 
derivedEJ. These excitations may be linked with hot spots created in peripheral 
heavy ion coUisionstJ. 

Let us consider the next type of nonlinear excitations in the region of nuclear 
surface, as a quadrupole deformation of the nuclear potential: 

U = UoScch\{r-R)/a)^2aUo^^^^^^^-^^^ (6) 

dr 

The evolutionjif this perturbation having the soliton form with time is given by the 
KdV equationc3: 

Ut - 6UU^ + PU^^^ = 0, (7) 
where the function U is the potential of Schroedinger equation 

13^^., -U^l> = 0, 13^ h^/2m (8) 

On the other hand this perturbation in the form of soliton can decay on the certain 
quantity of solitons with the small a mplitudcBil. If we fix the KdVE in the form 

Ut + UU^ + /JC/xxx = 0, (9) 

its solution will be in the form: 

U{x) = Uo sch^iUa/Upy^^X) (10) 
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Now we apply the similarity principle. The initial condition we take in the form: 

U{x,0)^Uo(l)ix,l). (11) 

With new variables 

V = U/Uo, x = X/l, T^Uot/l, (12) 
we obtain the new KdVE: 

Vt + mx + cr^^Vxxx = 0, (13) 

where a = 1{U(J3)^/'^. 

In the workEil the critical value a = \/l2 was obtained. With increasing number 
cr the corresponding perturbation breaks up into a larger number of solitons. The 
number of this solitons was determined by KarpmanE3 

iV = 2«^. (14) 

For the scattering of a-particles with kinetic energies in the range 4 and 12 MeV 
we can get for the nucleus ^^C the value N^5. It corresponds to the number of 
existing states generated by the surface potential well: K"^ — 0,2,4,6,8. Gardner 
et alr3 showed, that the amplitude of solitons determined by eigenvalues: v4 = 2i?„. 

The potential having the soliton shape gives the following spectrumcJ: 



En 



8mA 



(-(l-2n) + (l + 8mAt/o?i"^)^/^) . (15) 



For our example {n = 0,2,4,6,8) we can get the following values {MeV): —15.57, 
—8.65, —1.75, 5.19, 12.11 respectively. The smaller eigen value the lager soliton 
amplitude. Knowing the eigen values one can construct the total potential 0. 

This expression (6) provides the key to the understanding of the behaviour of 
the deformation lengthJor different targets and for different projectiles. In a recent 
paper by Koster et. al.E3 the empirical formula 

Z N 

{02fo)nucleo7is ~ ~^i(^2l"o) protons ~^ ~^i/^^''^o) neutrons i^^) 

is givfin to explain the different values of the deformation length for even isotopes 
of fM. 

In terms of the soliton approa,ch the compatibility condition for solitons in dif- 
ferent media can be writtenES 

A6 = A' 6' = const (17) 
where A and Ar correspond to the formula (17) and 

2/x 2fj.' 
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and 

The product AS should be the same for different isotopes of the target nucleus. 
The resulting prediction is that the deformation length in the case of inelastic 
scattering of a-particles on different Pd isotopes is should be largest for the heavest 
one for the 2+ state. 

Thus we interpreted two equations NoSE and KdVE respectively to the problem 
of clustering. This work is only the first approximation to this problem because we 
used only one dimension. To our mind the nextstep should be the application of 
Kadomtsev-Petviashvily equation to this problem. 

The most promising thing is the application of the Similarity principle to the 
linear alpha-chain states in nucleic^. 

4. Solitons as the clusters 

What are clusters? Clusters are spatially correlated nucleons. What kind of 
forces are between clusters and clusters and nuclei? Japanes theoreticians showedE3, 
that there are four typical " di-molecules" in nature (moleculas consisting from 
atoms, nuclei, nuclei plus hyperon, nucleons) with the same kind of potential. We 
consider the typical potential interaction between *He atoms. 

In the approximation one can consider He atoms as hard spheres with the radius 
2.7A. Minimum says that atoms have a tendency to develop into crystall. But the 
warm energy does not permit it. 

In the different energy scales the situations is very similar. 

The origination of the attractive parLpf the potential is due to the exchange of 
nucleons in the case of nuclear potentialEl. 

We shall consider the repulsive core. The best understanding can be obtained in 
the micro and macro approaches. The existence of a repulsive core in the potential 
interaction of two complex nuclei follows from Pauli principle which forbidds their 
mutual interpenetration. Let us give one example - the interaction of two alpha- 
particles. In this interaction the repulsive core is manifested as a result of which 
quasimolecular state the nucleus formes. The a + a phase shifts were obtained 
from the different experiments. The phase shift 62 begins to increase rapidly from 
an energy about 1 MeV and S4 from an energy about 6 MeV. This indicates 
an extension with a fairly abrupt boundary in the exterrior region with r < Afm 
(the impact parameter b = XI is approximately 4 fm at Ecm = 1 MeV and 3 fm 
at Ecm = 6 MeV ). This character of the scattering can be well described by a 
phenomenological potential with the repulsive core. 

In terms of the NoSE we can estimate the value of a hard core. For this purpose 
it is necessary to rewrite the NoSE. 



-— -r^U +{E- V)U + PU^ = 0. 
2m dx"^ 



(18) 
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After the factorization at the wave funtion U 

/ = 0, at X = 
/ = 1, at X ^ CO 



we obtain the NoSE 



^ + aij + (3tp^ ^ 0. (19) 

Zm ax'^ 



Then using such approaches a— >oo,[/ = 'i/;o,'0o = — a//? we have 

^^ + V = 0. (20) 



2ma dx^ a 

In the compact form this equation can be transformed into another equation 

Multiplaying this equatin on the derivative df / dx we can obtain a solution 

f ^th[x/\/2£]. (22) 

The solution is very similar on the correlation function of Boze gas with hard 
spheres: 

/ = 0, for r < 6, 
f^th[r/5~l]. 

Where 5 - is the radius of a hard sphere. In terms of this approach the radius of 
the hard core is: 

= V2^ = \/2— ^^^^^^^^=. (23) 

^2m{E - V) ^m{E - V) 

The microscopic estimates of the value of the hard core were made in the 
work of V.N.Bragin and R.DonangeloU and in the work of A.Tohsaki-Suzuki and 
K.NaitoLJ . The last athors used RGM. In this case it is very evident to expect 
the appearance of the hard core. In the first work the athors used the Brukner 
theory according to whitch the total energy of a system of interacting fermions can 
be written as a functional of the local single-particle energy density. They have 
computed the heavy ion potential in the sudden approximation as 

V{t) = J K(pi + P2) - e(pi) - C(P2)]dr . (24) 

The potential which was obtained in this approuch is highly repulsive in the 
nuclear interior quite an opposite behavior to that of the folding model calculations. 
Summing up, it is necessary to say, that the momentum and energy dependent hard 
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or soft core appears in the microseoptic calculations of two interacting nuclei due 
to Pauli prinsiple. On this base this core can be taken in account also in the 
macroscopic approaches. 

Keeping in mind the structure of the potential for nucleous-nucleous collision 
(the repulsive and attractive parts) we applied the inverse field method, which will 
be described in the next section. Within the inverse mean field method solitons 
are taken to model of elastic a + a collisions in a TDHF-like fasion. Attention is 
drawn to common points of this approach with TDHF. The analytical formula ijr 
the phase-shift within this approach yields a nice correspondence to experimentcil. 

5. Nuclear instability and soliton theory 

An existence of solitary waves is determined by two essential factors, namely, 
nonlinearity and dispersion. Both the factors, which are responsible for the stability 
of a wave, are connected in their turn with two different types of instability. A 
localized pulse will tend to spread out due to dispersion terms of the equations of 
motion. The nonlinearity which is responsible for the formation of solitary waves, 
on the other hand, automatically leads to their destruction, if it is alone. Both 
instabitities may compensate each other and lead to stable solutions (solitons) . 

Let us look from these points of view at multifragmentation phenomena. The for- 
mation and breakup of a highly excited and compressed nuclear system is the most 
striking process observed in intermediate-energy heavy ion reactions. How does such 
a system expand and finally disassemble when passing through a regime of dynami- 
cal instabilities? What is the mechanism of the clustering (stable light and interme- 
diate mass fragment production)? Quite a variety of models have been developed to 
discuss this question (see, for instance, the recent reviewEa and ProceedingstZl) . The 
dynamical clusterization in the presence of instabilities is the focus of attention of 
the intermediate energy heavy ions physicsPfr^. Ten years agojnultifragmentation 
has been associated with the onset of the spinodal instabilityc3. This instability 
is associated with the transit of a homogeneous fluid across a domain of negative 
pressure, which leads to its breaking up into droplets of denser liquid embedded in a 
lower density vapor. Since the spinodal instability can occur in an infinite system, it 
can be called the bulk or volume instability. On the other hand, it physically means 
that pressure depends on density, that is just a nonlinearity in terms of density H. 

Recently^, it has been pointed out that a new kind of instability (sheet instabil- 
ity) may play an important role in multifragmentation. This new instability can be 
assigned to the class of surface instabilities of the Rayleigh- Taylor kindEl. System 
escapes from the high surface energy of the intermediate complex by breaking up 
into a number of spherical fragments with less overall surface. At the same time, it 
physically means the existence of the gradient terms of the equations of motion, i.e 
the dispersioncll. 

The spinodal instability and the Rayleigh- Taylor instability may compensate 
each other and lead to stable quasi-soliton type objects. In the next Sec. we 
present this physical picture using a simple analytical model proposed to describe 
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the time evolution of compressed nuclear systems. 

An analysis of stability of nonlinear dynamical systems and an analysis of non- 
linear evolution of initial complex states is a traditional goal of Soliton Theory. The 
inverse methods to integrate nonlinear evolution equations arc often more effective 
than a direct numerical integration. Let us demonstrate this statement for a very 
simple casecj. 

The type of systems under our consideration are uncharged slabs of nuclear 
matter. The slabs are finite in the z coordinate and infinite and homogeneous in 
two transverse directions. 

The basic equations for the static mean-field description of the slabs are the 
following 

1 h'^P 

i/'kin(x) = — =V'«(2)exp(ikj_r), ekj_„ = -z-^ + e„, 

No 

p{x) =^ p{z) = an^pl{z), (25) 

n=l 

No r, 

2m 



E 
A 



A =4> A = (6A/9^/7r)^/^ = a„, a„ = ^^(eF - e„), 

n— 1 

^0 ^oo J,/. _ No 



[p{z)]dz, 



n=l n=l 

where r = {x,y),hj_ = {kx,ky), and Q is the transverse normalization area. A^o is 
the number of occupied bound orbitals. 

The dynamical description will be done in the framework of the inverse mean 
field method. One can found the details of this approach i We concentrate 

here only on essentials. 

The evolution of a system is given by the famous hydrodynamical Korteweg-de 
Vries equation (KdV) for the mean-field potential U{z,t) 

dU _ ^^^dU d^U 

n — l ^ ' 

where Sn are constants which are determined by the initial conditions. 

General solution of KdV Eq. (26) can be derived in principle via direct methods 
numerically. This way is to assign a functional of interaction £ (as usual an effective 
density dependent Skyrme force), a total number of particles (or a thickness of a 
slab A) and to solve Hartree-Fock equations to derive a spectrum of the single 
particle states e„ and wave functions ■0„(z,O), the density profile p{z,0) and the 
one-body potential U{z,0) for the initial compressed nucleus. Then, one calculate 
an evolution of the one-body potential with the help of Eq. (26). 
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However, there is an inverse method to solve KdV Eq. (26). The main advantage 
of this way is to reduce the solution of the nonlinear KdV Eq. (26) to the solution 
of the linear Schroedinger - type equation 

- — --j^i:n{z,t) + U{z,t)ilj„{z,t) = e„ilj„{z,t), (27) 

and linear integral Gelfand - Levitan - Marchenko equation to derive the function 
K{x,y) 

K{x,y)+B{x + y)+ / B{y + z)K{x, z)dz = 0. (28) 

The kernel B is determined by the reflection coefficients R{k){ek = h^k"^ /2m). and 
by the N bound state eigenvalues e„ = — ?i,^K^/2m. N is the total number of the 
bound orbitals. 

B{z) = y C2(k„) + - / R{k)exp{ikz)dk. 
The coefficients C„ are uniquely specified by the boundary conditions 

Cn(Kn) = lim 1pn{z)exp{KnZ), 

z—*oo 

and the wanted single particle potential is given by 

U{z,t) = --^K{z,z). 
m oz 

The time t is included in Eqs. (27,28) only as a parameter, so it has been omitted 
in the above formulas. 

The general solution, U{z,t), should naturally contain both, contributions due 
to the continuum of the spectrum and to its discrete part. There is no way to obtain 
the general solution U {z, t) in a closed form. Eqs. (27,28) have to be solved only 
numerically. 

However in the case of reflectless {R{k) = 0) symmetrical {U{—z,0) = U{z,0)) 
potentials one can derive the following basic relations 

^(^'*) = ^Mdet ||M||) = Yl ^n€{z,t), 

TDj (jZ Tft 

n=l 

N 



tPn{z,t) = Y,iM-%Mz,t), 



n=l 



Xn{z, t) = Cn{Kn) exp{-KnZ + 2h'^ K^Snt/m)), (29) 
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Mni(Z,t) = dnl H 



N 

Cn{Hn) — I 2Kn I I 



1/2 



So, the wave functions, potential and the density profile are completely defined by 
the bound state eigenvalues. 

The first step is to solve the Schroedinger Eq. (27) for the initial potential 
U{z,Q), which is suitable to simulate compressed nuclear system or to simulate this 
state with the help of spectrum. Then one calculates the evolution of p{z,t) and 
U{z,t) with the help of Eqs. (29). 

Although there is definitely some progress in the application of the inverse meth- 
ods to nuclear physics they are not yet too popular. As illustration of these methods, 
we considered a one-dimensional three-level system in detailsEl. 

A three- level system may be useful for modelling the evolution of light nuclei, 
for instance, of oxygerH. 

for large z and t the time - dependent one - body potential and the corresponding 
density distributions are represented by a set of stable solitary waves. The energy 
spectrum of an initially compressed system completely determines widths, velocities 
and the phase shifts of the solitons. 

The number of waves is equal to the number of occupied bound orbitals. Thick- 
ness ('number' of particles) of an n - wave is equal to a„. 

Reflecting terms (i?(fc) 7^ 0) of U{z,t) cause ripples (oscillating waves of a small 
amplitude) in addition to the solitons in the final state. 

The initially compressed system expands so that for large times one can observe 
separate density solitons and ripples ('emissions'). This picture is in accordance with 
the TDHF simulation of the time evolution of a compressed O^® nucleustll. The 
disassembly shows collective flow and clusterization. It is important to note that 
the clusterization was not observed in the absence of the self-consistent mean-fleld 
potential, i.e this confirmes our supposition that the nonlinearity is very important 
for the clustering. 

It is necessary to note that the present model is too primitive in order to describe 
a real breakup process. However this model can be used to illustrate an inverse 
mean-field method scheme, a nonlinear principle of superposition and the idea that 
nonlinearity and dispersion terms of the evolution equation can lead to clusterization 
in the final channelt2l. 

Certainly there are a lot of open questions relative to the presented approach. 
The most crucial concerns the generalization to 3-1-1 dimensional model with a finite 
temperature. One possible way to do it keeping symplicity of the ap pro a.ch would 
be to use an information theoryS. Such investigations are in progress 
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6. Summary 

It is shown that very different effects of nuclear cluster physics could be explained 
in the frame of Korteveg de Vries Equation and cubic Nonlinear Schrodinger Equa- 
tion. 

The NoSE has been shown to be a very useful tool for the analysis of experimen- 
tal data involving the clastic interaction of light and heavy ions. Its novel future is 
that it determines fairly accurately the compressibility modulus of nuclear matter. 
The problem of dynamical instability and clustering are considered from the points 
of view of the soliton conception and using KdV type equation for a mean-field 
potential. Nonlinear nuclear vortex is discussed.. Due to the soliton conception it is 
understandable the behavior of the quadrupole deformation parameters for different 
isotopes of one nucleus. 
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